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ABSTRACT 

Observational measurements of active region emission measures contain clues to 
the time-dependence of the underlying heating mechanism. A strongly non-linear scal- 
ing of the emission measure with temperature indicates a large amount of hot plasma 
relative to warm plasma. A weakly non-linear (or linear) scaling of the emission mea- 
sure indicates a relatively large amount of warm plasma, suggesting that the hot active 
region plasma is allowed to cool and so the heating is impulsive with a long repeat 
time. This case is called low-frequency nanoflare heating and we investigate its fea- 
sibility as an active region heating scenario here. We explore a parameter space of 
heating and coronal loop properties with a hydrodynamic model. For each model run, 
we calculate the slope a of the emission measure distribution EM(T) oc T a . Our 
conclusions are: (1) low-frequency nanoflare heating is consistent with about 36% of 
observed active region cores when uncertainties in the atomic data are not accounted 
for; (2) proper consideration of uncertainties yields a range in which as many as 77% 
of observed active regions are consistent with low-frequency nanoflare heating and 
as few as zero; (3) low-frequency nanoflare heating cannot explain observed slopes 
greater than 3; (4) the upper limit to the volumetric energy release is in the region of 
50 erg cnT 3 to avoid unphysical magnetic field strengths; (5) the heating timescale 
may be short for loops of total length less than 40 Mm to be consistent with the ob- 
served range of slopes; (6) predicted slopes are consistently steeper for longer loops. 

Subject headings: Sun: corona 
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1. Introduction 

One of the most enduring problems in modern astrophysics is to explain how the million- 
degree solar corona is created and sustained. The corona is a highly non-uniform environment 
which reflects its magnetic structure and activity. For example, there are magnetically open and 
closed regions; in open regions the field lines extend great distances into the heliosphere where 
they may reconnect with planetary magnetic fields, and in closed regions magnetic loops, line-tied 
at the solar surface and illuminated by plasma at EUV and X-ray temperatures, are observed. 
Quiet Sun regions are associated with relatively weak magnetic fields, whereas active regions, 
which overlie the sunspots at the surface, have stronger magnetic fields and give rise to a great 
deal of dynamic and eruptive solar activity. The magnetic fields are strongest in active region 
cores, where temperatures can exceed 5 MK. Regions of weaker magnetic field are generally 
associated with lower temperatures and it therefore seems clear that the heating rate must be 
in some way related to the magnetic field. The most obvious connection is through heating 
the plasma via the release of magnetic energy. This might occur by gradually stressing the 
magnetic fields threading the corona, building up the energy stored in the field, until the energy 
is released by reconnection leading to direct heating, bulk motion, and p article acceler ation, 



the latter two ultimately being thermalized via collisional processes (e.g. 



Parker 



1988), for 



example. Alternatively, Alfven waves propagating along t he field lines may inter a ct, leading t o 



energy release and dissipation via resonant processes (e.g. 



Mcintosh et al. 



2011b . 



Klimchuk 



(|2006r) provides a more detailed discussion of these and other possibilities. Unfortunately, direct 



observational signatures of co ronal heating are extremely difficult to detect (Bra dshaw & Cargill 



2006; 



Reale & Orlando 



2008) and so the actual mechanism which raises the coronal plasma to 



such tremendous temperatures has not yet been identified. 



Nonetheless there are observational measurements that may provide information concerning 
the properties of the heating mechanism, such as the timescale over which the release and 
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dissipation of energy must occur, which may in turn yield clues to the physical nature of the 
mechanism itself. One popular conception of coronal heating is to consider the manner in which 
a monolithic (as observed) coronal loop is heated; if one assumes that the loop is composed of 
many thermally isolated, sub-resolution filaments, then the question to be answered concerns how 
each individual filament is heated. A single filament is understood to be the thinnest magnetic flux 
tube with an isothermal cross-section. There are two broad possibilities for the time-dependence 
of the heating mechanism for each filament; either steady or impulsive. In the case of steady 
heating energy is released and dissipated at a more or less constant rate, and the mechanism 
operates for a period of time that is much longer than the cooling timescale. In the case of 
impulsive heating energy is released and dissipated on a timescale that is significantly shorter 
than the cooling timescale. In the former case the filament heats up and eventually reaches some 
new, hot, hydrostatic equilibrium. In the latter case the filament heats up and is then allowed to 
cool and drain. We will refer to impulsive heating events by the generic term nanoflare, which 
is understood to mean any mechanism that gives rise to the impulsive release and dissipation of 
energy. 

An important parameter is the frequency with which nanoflares occur on the same filament. 
For example, if a filament is heated and then cools and drains, returning to its initial state before 
being re-heated, then we refer to this as low-frequency nanoflare heating. However, if a filament 
is heated and then re-heated just a short while into its cooling and draining phase, then we refer 
to this as high-frequency nanoflare heating. We note that steady heating is just the upper limit 
of high-frequency nanoflare heating, where the delay between successive nanoflares on a single 
filament is effectively zero. 

There has been a great deal of recent interest in one potential diagnostic of the frequency 
of occurrence of heating events in active region cores: the gradient or slope (a) of the 
emission measure (in log - log space) between the temperature of peak emission measure and 
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1 MK ( 


Warren et a 




201 


1 


Schmelz & Pathak 


2012 


). 



it IWinebargeretal. 11201 ll : iTripathi et al. 11201 ll : IWarren et al. 112012 : 



2012|) . Of the active region cores analysed to date the observationally 
measured slopes are in the range 1.70 < a < 5.17. Steeper slopes indicate a greater proportion 
of hot (e.g. 3-5 MK) material relative to warm (« 1 MK) material and this is consistent with 
steady heating (truly steady or high-frequency nanoflare heating). In this scenario most of the 
sub-resolution filaments are maintained at high temperatures, with relatively few cooling and 
draining; consequently there is little warm material. Shallower slopes indicate commensurately 
more equal proportions of hot and warm material, consistent with low-frequency heating. In this 
scenario the filaments are allowed to cool, with the result that there is more warm material than in 
the steady heating limit. 

Recent work has focused on determining whether low-frequency nanoflares are consistent 
with the observed emission measure slopes and what the properties of the individual nanoflares 
(e.g. volumetric heating rate, duration) heating each fila ment must be, as we ll as the properties 



Warren et al. 



(120111) found that neither 



of the loops themselves (e.g. initial conditions, length), 
low-frequency nanoflares nor steady heating could explain their observed emission measure 
slopes (« 3.26), tho ugh later observations yielded significantly shallower slopes (2.05 - 2.70: 



Tripathi et al. 



20111) which were more consistent with low-frequency heating. Steady heating 
led to effectively isothermal emission measures which severely underestimated the amount 
of warm plasma pr e sent, especially in diffuse regions where no discrete loops are observed. 



Mulu-Moore et al. 



(2011) carried out several numerical experiments to study the emission 
measure slope in the case of heating by low-frequency nanoflares, focusing in particular on the 
influence of the radiative losses by using different abundance sets in their models. They found 
slopes in the range 1.60 < a < 2.00 for photospheric abundances and 2.00 < a < 2.30 for a 
coronal abu ndance set with a n enh anced population of low first ionization potential (FIP) elements 



(e.g. iron). 



Reale & Landi 



(|2012r) have also discussed the influence of radiative losses on the 



emission measure slope; using the most recent atomic data they found steeper emission measure 
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slopes below 2 MK, which they associated with enhanced cooling in the 1-2 MK range. 

In the present work we focus on the heating of active region cores by low-frequency 
nanoflares and carry out an extensive survey of the parameter space of possible heating events 
and loop properties. In particular, we consider the importance to the emission measure slope of 
the volumetric heating rate, the duration of the heating t h , the temporal envelope of the heating 
profile, the initial conditions in the loop (n, T), and the loop length. We carry out extremely 
detailed forward modeling to produce synthetic spectra and use this to calculate the emission 
measure for each of our numerical experiments. In Section [2] we describe our modeling approach 
and the experiments there were carried out. We discuss our results in Section [3] and, finally, present 
a summary of our key results and a number of conclusions concerning heating by low-frequency 
nanoflares in Section HI 



2. Numerical model and experiments 



Our approach to carryin g out the numerical experimen ts and our forward modeling procedure 



are documented in detail in 



Bradshaw & Klimchuk 



(|201l|) . However, we will provide a brief 



summary and then focus on what is new here. We have explored an extensive parameter space 
of coronal loop and low-frequency nanoflare properties that we believe may be representative 
of active region core heating, and we present the results from 45 model runs. A low-frequency 
nanoflare is defined such that: (a) the heating timescale is significantly shorter than the timescale 
required for the strand to reach a new hydrostatic equilibrium if the increased heat input had 
remained steady; and (b) the time between each individual heating event on a single strand must 
be significantly longer than the cooling timescale following heating. 

The first five columns of Table Q] summarize the key loop properties and heating parameter 
values for each run. The first column indicates the number assigned to each Run, which will 
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be used to reference them as they are discussed in turn. The second column gives the loop 
length 2L defined as the distance between the transition region foot-points along the loop. The 
computational domain also includes a deep (many scale heights) chromosphere at 2 x 10 4 K 
that is attached to each transition region foot-point in order to provide a source of materia l to 



ablate into the corona upon heating and to maintain a stable atmosphere (jPeres et al 



1982). The 



third and fourth columns describe the properties of the impulsive heating event. E H0 and t h are 
the peak volumetric heating rate and the duration of heating. The fifth column gives the total 
volumetric heating input E H . Note that the temporal envelope of the heating profile can be square 
(E H = E H0 x t h ) or triangular (E H = 0.5 x E m x r H ). The first 14 entries in Tabled] correspond 



to the 14 numerical experiments that form the basis of 



Bradshaw & Klimchuk 



(120111) and the 



remainder were conducted to substantially broaden the parameter space. 

Table 1:: A Summary of the Numerical Experiments Relat- 
ing to the Loop and Impulsive Heating Event Properties. 



Run# 2L E m t h E h log 10 T peak a mo d e l ^observed 





[Mm] 


[erg cm 3 s '] 


[s] 


[erg cm 3 ] 


[K] 






1 


20 


0.05 


10 


0.5 


5.85 






2 


20 


0.10 


10 


1.0 


5.95 






3 


20 


0.10 


30 


3.0 


6.15 


0.58 


1.17 


4 


20 


0.10 


100 


10 


6.45 


0.76 


0.81 


5 


20 


0.10 


300 


30 


6.35 


0.98 


0.83 


6 


20 


1.00 


10 


10 


6.35 


1.38 


0.89 


7 


20 


1.00 


30 


30 


6.45 


1.95 


1.65 


8 


20 


5.00 


10 


50 


6.55 


1.79 


1.73 


9 


80 


0.10 


10 


1.0 


6.25 


1.05 


1.69 


10 


80 


0.10 


30 


3.0 


6.45 


1.21 


1.29 
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Table 1:: A Summary of the Numerical Experiments Relat- 
ing to the Loop and Impulsive Heating Event Properties. 



Run # 


9f 




T TT 




10 gio J peak 


"model 


"observed 




rMml 

|^±VA111J 


L^lti V^lll l5 i 




L cl fe ^ m J 


L^J 






1 1 


80 

O V 


10 


100 

1 \J\J 


10 


6 45 


9 94 


9 99 


1 9 


OV 


10 

u. 1 u 


300 

juu 


30 

JU 


6 6S 


1 Q0 


1 8S 

1 . J 


1 3 


80 


1 00 
i .uu 


10 

1U 


10 

1U 


6 SS 

U.J J 


1 8Q 

1 .0:7 


9 19 

Z. 1 z 


1 A 

It- 


80 
ou 


1 00 


30 

JU 


30 

JU 


6 6S 


9 01 

Z.Ul 


1 80 
l .OU 


1 < 

1 J 


AO 


O 03 
U.UJ 


soo 

JUU 


7 S 
/ .J 


6 AS 
U.'+J 


1 1 A 
1 . it 


1 97 
l.Z / 


1 6 
1U 


AO 

H-U 


06 

U.UU 


SOO 

JUU 


1 S 

l J 


6 SS 

U.J J 


1 39 

1 . JZ 


1 68 
l .UO 


1 7 


40 

M-U 


O 10 
u. iu 


SOO 

JUU 


9S 

Z J 


6 SS 

U.J J 


1 S8 
1 .Jo 


1 1 S 

1 . 1 J 


1 8 


AO 


SO 


soo 

JUU 


1 9S 

1ZJ 


6 7S 
U. / J 


1 QS 

1 .yJ 


1 03 
1 .y J 


1 Q 

Ly 


80 


03 


soo 

JUU 


7 S 
/ .J 


6 SS 

U.J J 


1 AO 


1 8S 
1 .0 J 


90 

ZU 


80 


06 

U.UU 


SOO 

JUU 


1 < 

1 J 


6 6S 

U.UJ 


1 S3 
1 .J J 


9 1 6 

Z. 1U 


91 
Zl 


80 

ou 


10 
U. lu 


SOO 

JUU 


9S 

ZJ 


6 6S 

U.UJ 


1 73 


1 SA 


99 

ZZ 


80 

ou 


SO 

U.JU 


SOO 

JUU 


1 9S 

1ZJ 


6 QS 

0.7J 


1 79 
1. /z 


9 OQ 
z.uy 


93 


1 60 
1UU 


03 
U.UJ 


soo 

JUU 


7 S 
/ .J 


6 6S 
U.UJ 


1 68 
1 .Do 


9 A7 


24 


160 


0.06 


500 


15 


6.65 


2.01 


1.77 


25 


160 


0.10 


500 


25 


6.75 


1.90 


1.95 


26 


160 


0.50 


500 


125 


7.05 


1.74 


2.20 


27 


160 


0.03 


2000 


30 


6.75 


1.97 


2.00 


28 


160 


0.06 


2000 


60 


6.85 


1.92 


2.14 


29 


160 


0.10 


2000 


100 


7.15 


1.62 


2.23 


30 


160 


0.50 


2000 


500 


7.35 


1.71 


2.26 


31 


40 


0.03 


500 


15 


6.45 


1.45 


1.66 
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Table 1:: A Summary of the Numerical Experiments Relat- 
ing to the Loop and Impulsive Heating Event Properties. 



Run # 

1VU11 IT 






T TT 




lu gio J peak 


"model 


"observed 




rMml 


trover ptn - ^ c ^1 




L ei g Lm 1 








19 


40 


06 
U.UU 


500 


10 


6 65 

U.UJ 


1 49 


1 94 


11 


40 


10 


500 

J uu 


50 

JU 


6 65 


1 61 

1 .U J 


1 50 

1 . JU 


14 

JT" 


80 


01 

U.UJ 


500 

J uu 


i j 


6 55 

U.J J 


1 81 


9 16 


15 


80 
ou 


06 

U.UU 


500 

JUU 


10 


6 65 

U.UJ 


1 78 
l . / o 


1 .uu 


16 
JU 


80 

OU 


10 

U. 1 U 


500 

JUU 


50 

JU 


6 75 

U. / J 


1 81 

1 . 1 


1 R? 


j i 




0^ 

U.UJ 


son 


1 J 


6 6S 


1 QR 


1 . / u 


38 


160 


0.06 


500 


30 


6.75 


2.03 


2.04 


39 


160 


0.10 


500 


50 


6.85 


1.93 


2.14 


40 


80 


0.20 


10 


1.0 


6.35 


1.49 


2.64 


41 


80 


0.20 


30 


3.0 


6.45 


1.46 


2.56 


42 


80 


0.20 


60 


6.0 


6.45 


1.91 


1.93 


43 


80 


0.20 


120 


12 


6.55 


1.95 


2.28 


44 


80 


0.20 


300 


30 


6.65 


2.05 


1.90 


45 


80 


0.20 


600 


60 


6.75 


2.04 


2.16 



We solve the one dimensional hydrodynamic equations appropriate for a single magnetic 
strand in the field-aligned direction. The diameter of the strand cross-section is assumed to be so 
small that many strands are visible along the line of sight for any observing instrument. In the 
case of a discernable loop, the many strands must combine in the resolved volume to give rise to 
the observed loop structure. We adopt a multi-species approach by treating electrons and ions 
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as separate fluids, and couple them via Co ulomb collisions. We preferentially heat the electrons 



(though see 



Longcope & Bradshaw 



2010, for an alternative scenario) and therefore expect that 
T e ± T t during the heating and conductive cooling phases, and that T e « T, at the onset of radiative 
cooling. We assume quasi-neutrality (n e = n { = n) and current free (v e = v, = v) conditions. 

The equations solved are formulated to describe the conservation of mass, momentum and 
energy. They include transport and compression, viscous stress, gravitational acceleration and 
potential energy, Coulomb collisions, thermal conduction, optically-thin radiation, and external 
energy input (heating). Collisions between like- species are frequent enough that the electron and 
ion equations of state are given by p e = k B nT e and p, = k B nTj, and only the thermal component 
of the electron energy and the thermal plus kinetic components of the ion energy are considered. 
Viscous interactions are expected to become extremely important at the hig h temperatures r eached 



by some of our experiments and so must be included in the strand physics (jPeres & Reale 



1993) 



The optically-thin radiation calcul ation accounts for the volumetric energy loss due both to line 



emis sion (for corona l abun dances 



Feldman et al. 



v6 



Dere et al. 



1991 



19921) and thermal bremsstrahlung (Chianti 



2009). Runs 1-14 and 40 - 45 account for the effect of non-equilibrium 
ionization of Ca and Fe on the radiative losses, but deviations from equilibrium are only significant 
at times when radiation plays a very minor role in the energy balance (e.g. during heating and 
conductive cooling) and so the effect on the plasma thermodynamics is small. Nonetheless, 



the effect on the emission spectrum can be significant at hig 



temperature of peak emission measure: 



ler temperatures (e.g. above th e 



Bradshaw & Cargill 



2006 



Reale & Orlando 



2008) 



The spatial profile of the volumetric energy input (heating) is such that it is uniform along the 
magnetic field and we adopt one of two temporal envelopes: (a) constant heating for a period t h 
(Runs 1-14 and 31 - 39); and (b) a triangular profile with a linear increase to the maximum 
heating rate at r H /2 and then a linear decrease to zero at r H (Runs 15-30 and 40 - 45). 



We consider only spatially uniform heating in our numerical experiments for two reasons: 
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(a) we are wary of introducing an additional, unconstrained parameter into our study; and (b) 
we believe that the radiative and enthalpy-driven cooling and draining phase, which determines 
the emission measure, is relatively insensitive to the details of the spatial distribution of heat 
in the impulsive case. Our justification for (b) is based on the efficiency of thermal conduction 
at redistributing the energy released into the loop. In order to obtain emission measures that 
peak at temperatures (T^^) that are consistent with the observed range, it is necessary to heat 
the strand to temperatures significantly greater than T^ ea ^ (e.g. into the region of 10 MK). At 
these temperatures, and at the rarefied densities of the initial conditions, thermal conduction 
redistributes the energy throughout the loop in just a few seconds, which is effectively akin to 
uniform heating on the timescales of interest to us. By the time the strand enters the radiative 



cooling a nd draining phase then all 'memory' of the initial spatia 



been lost (IWinebarger & Warren 



2004; 



Patsourakos & Klimchuk 



distribution of heat will have 



2005) 



We use the numerical code HYDRAD (IBradshaw & Mason 



2003 



Bradshaw & Klimchuk 



201 1|) to solve the two-fluid hydrodynamic equations under the circumstances described above. 
HYDRAD has several desirable features which make it ideal for application to the study of 
impulsive heating to extremely high temperatures. Written exclusively in C++, it is fast and 
robust, and models an entire loop strand (foot-point to foot-point for any geometry via an analytical 
equation, or look-up table, for gravity) with an adaptive grid for efficiently capturing small-scale 
properties of the solution. It is user-friendly and easily configurable via a Java-developed graphical 
user interface. 



The forward mode l ing a spect of this work follows the procedure described by 



Bradshaw & Klimchuk 



(|201l|) . Although the EUV emission from active region cores is 



generall y diffuse (only a fractio n of the EUV emission is contained in observationally discernable 



loops: IViall & Klimchuk 



201 lh . we model the core as though it were a single multi-stranded 



loop. We constructed a snapshot of the loop from the data output by each experiment. The loop is 
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comprised of sub-resolution magnetic strands, where each strand represents one stage (captured in 
1 second intervals) of the heating and cooling cycle. Therefore, a loop with a heating and cooling 



Guarrasi et al. 



2010) 



cycle lasting 2000 s would be composed of 2000 individual strands (see also 
We then calculated the emission measure in the region of the loop apex, since we are interested in 
active region core heating and wish to avoid foot-point / moss contamination, using two different 
methods. The first gives the 'true' or model emission measure, which is obtained directly from the 
numerical data (n and T) stored in each grid cell. Th e second is the ' observed' emission m easure, 



forward mod eled by applying th e Pottasch method (Pottasch 



1963 



Hinode-ElS (Cu lhane et al. 



Tripathi et al. 



2011) to a 



|2007|) spectrum synthesized from the numerical data. The intensity 
of emission from each resolved volume element along the loop, recorded on the corresponding 
detector pixel, is calculated by summing the line-of-sight emission from each sub-resolution 
strand fold ed through the appropriate inst rument response function (e.g. Equations (8), (9) and 



Figure 1 of iBradshaw & Klimchuk 



201 1|) . Note that the 'observed' emission measure is idealized 



in that there are no errors due to photon statistics, photometric calibration of the instrument, 
or incorrect atomic physics (abundance, collision and excitation rates, etc.). The same atomic 
physics is used to generate the synthetic spectrum as is use d to infer the emission measures from 



that spectrum, so any errors cancel out. 



Landi et al. 



(120121) have investigated the effect of using 



inconsistent atomic physics in differential emission measure (DEM) diagnostics and found that 
it can lead to uncertainty in the peak DEM value. Uncertainties associated with real data are 
substantial and can have a significant impact on the emission measure slope, as we discuss later. 

We performed this calculation for a number of individual spectral lines and applied the 
Pottasch method to each (as though they were observed lines) in order to derive an emission 
measure plot. An example is shown in Figured] for Run 44. The plus signs show the observed 
measures and the diamonds show the true values. EM loci curves are also plotted for each line. 
These curves would intersect at a single temperature if the plasma were isothermal, and the 
degree of deviation from a single intersection point is an indication of the multi-thermality of the 
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plasma (ILandi & Klimchuk 11201 0|) . We used a linear least-squares fit to calculate the slope of each 
emission measure distribution between the temperature of its peak and 10 6 K. A steeper slope 
is consistent with a narrower, more isothermal, emission measure and a shallower slope with a 
broader, multi-thermal distribution of plasma. The spectral lines in the EIS wavelength channels 
used to construct the forward modeled emission measures used for our study are listed in Tabled 



Table 2:: A List of the Spectral Lines Used by the For- 
ward Modeling Code to Generate the Emission Measure 
Loci Plots. 



Ion Wavelength (A) log 10 T (K) 



Mg V 


i /o.j /y 


C A C 

J.4D 


Mg VI 


268.991 


5.65 


Mg VI 


270.391 


5.65 


Si VII 


275.354 


5.80 


Mg VII 


278.404 


5.80 


Mg VII 


280.745 


5.80 


FelX 


188.497 


5.85 


FelX 


197.865 


5.85 


Si IX 


258.082 


6.05 


FeX 


184.357 


6.05 


FeXI 


180.408 


6.15 


FeXI 


188.232 


6.15 


SiX 


258.371 


6.15 


SiX 


261.044 


6.15 


sx 


264.231 


6.15 


FeXII 


192.394 


6.20 
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Table 2:: A List of the Spectral Lines Used by the For- 
ward Modeling Code to Generate the Emission Measure 
Loci Plots. 



Ion 


Wavelength (A) 


log 1fl T (K) 

olU V / 


FeXII 


195.119 


6.20 


FeXIII 


202.044 


6.25 


FeXIII 


203.828 


6.25 


FeXIV 


264.790 


6.30 


FeXIV 


270.522 


6.30 


FeXIV 


274.204 


6.30 


FeXV 


284.163 


6.35 


SXIII 


256.685 


6.40 


FeXVI 


262.976 


6.45 


Ca XIV 


193.866 


6.55 


Ca XV 


200.972 


6.65 


CaXVI 


208.604 


6.70 


Ca XVII 


192.853 


6.75 


Fe XVII 


269.494 


6.75 



3. Results 

The key results of our study are presented in the final three columns of Tabled] The sixth 
column of Table Q] gives the temperature of the emission measure peak (log 10 7^^) and the 
final two columns give the slopes (a) of the model and forward modeled emission measures for 
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6.0 < log 10 T < log 10 7"p ea k, where EM(T) oc T a . The first finding to which we draw attention are 
the differences between the slopes calculated for the model and the forward modeled 'observed' 
emission measure. We find a general tendency for the slope in the forward modeled case to be 
steeper than the model slope as the loop length increases; it is true for 33% of the 20 Mm loops, 
43% of the 40 Mm loops, 68% of the 80 Mm loops, and 82% of the 160 Mm loops in Tabled] 
This implies that the slope calculated from observationally derived emission measures may not 
reflect the true slope, which would be the slope obtained if one were in possession of 'perfect' 
data (e.g. in the model case this would be the actual grid cell densities and temperatures). Forward 
modeling then allows us to investigate how the biases introduced by commonly used observational 
tools and techniques may influence sets of results and the conclusions associated with them. 



Warren et al. 



toil) 



By way of an example consider the Mg VII 280.745 A line, which 
noted was not consistent with other lines formed at a similar temperature (e.g. Si VII 275.354 A). 
They conjectured that since the Si X lines in their study were consistent with the Fe lines then the 
issue must lie with the Mg abundance and it was multiplied by a factor of 1.7 to force agreement 
between the Mg VII 280.745 A and the Si VII 275.354 A lines. However the upper panel of 
Figure \T\ shows an example emission measure loci plot, calculated using a spectrum forward 
modeled from numerical data, with two lines towards the low temperature range that are not 
consistent with their neighbours. The line formed at 10 5 8 K is the Mg VII line, but the line formed 
at 10 6 05 is the Si IX 258.082 A line. Since the Si IX line is not consistent with the neighbouring 
Si X lines then the problem cannot be related to the element abundance. Instead, it is related to the 
density sensitivity of these emission lines. The emission measure loci plot in the upper panel of 



Figure Q] was calcu 



(e.g. 



Tripathi et al. 



ated for a density of 10 cm , commonly assumed for an active region core 



l201l|) . in the contribution function. The lower panel of Figure Q] shows the 
same emission measure loci plot calculated for a density of 10 9 cm -3 in the contribution function. 
The Mg and Si lines are brought into consistency with one another and with the lines that are not 
density sensitive (which do not change). In addition, there is less scatter between lines formed at 
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Fig. 1. — Emission measure loci plots for Run 44. Diamonds are the 'perfect' model emission 
measure. + signs are the emission measure values calculated by applying the Pottasch method 
to a forward modeled Hinode-EIS spectrum. Top (bottom) plot assumes a density of 10 10 cm 3 
(10 9 cm -3 ) for the contribution functions. Pixel 42 is the loop apex pixel on the virtual detector. 
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the same temperature. The improved agreement between the density sensitive and non-density 
sensitive lines is a strong indication that the lower density is a better estimate of the active region 
core density. The best value is probably somewhat greater, but closer to 10 9 cm -3 . 



Tripathi et al. 



(2011) 



In the current work we use a value of 10 10 cm 3 for consistency with 
who adopt this value and also employ the Pottasch method for calculating their emission measures. 
Our aim is to be as consistent as possible with observational studies in our forward modeling. We 
note that using different densities in the contribution function will yield different emission measure 
slopes. The differences between the model and forward modeled slopes is of particular relevance 
to studies of emission measure slopes as diagnostics of the time-dependence of the underlying 
heating mechanism. Previous studies have found that emission measure slopes calculated from 
'perfect' model data are too shallow, in the case of low-frequency nanoflares, to explain observed 
slopes in excess of 3. 

Table 3:: A Summary of Recent Observational and Numeri- 
cal Modeling Results Concerning the Slope of the Emission 
Measure Coolward of the Peak. 



Slopes 



Range of log 10 T 



Reference 



Comments 



3.26 



2.17, oo 



6.00 - 6.60 



Warren et al. ("20111 



1 AR; 10" x 15" sub-region; 
MCMC method; with background 

Model slopes; 
low & high frequency nanoflares 



3.20 



6.00 - 6.50 



Winebarger et al. (201 1} 



1 AR; 5" x 25" sub-region; 
xrt_dem_interative2 .pro ; 
background subtract 



2.08 - 2.47 5.50 - 6.55 
2.05 - 2.70 



Tripathi et al. (20111 



with background 
background subtract 
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Table 3:: A Summary of Recent Observational and Numeri- 
cal Modeling Results Concerning the Slope of the Emission 
Measure Coolward of the Peak. 



Slopes 


Range of log 10 T 


Reference 


Comments 


2 ARs; 5" x 5" to 10" x 15" sub-regions 


1.60-2.00 


6.00 - [6.60 - 6.80] 


Mulu-Moore et al. (20 11) 


Model slopes; photospheric abundances 


2.00 - 2.30 






coronal abundances 


1.70-4.50 


6.00 - 6.60 


Warren et al. (2012) 


11 > 3.00,5 « 2.00; 



12 ARs; 2 with 2 sub-regions; 
1 with 3 sub-regions; multiple pixels; 
MCMC method; with background; 
slope increases with unsigned flux 



1.91 -5.17 6.00 - [6.30, 6.80] Schmelz & Pathak (2012) 4 < 2.60, 2 > 3.00; 

5 ARs; 2 with 2 sub-regions; 
3 with 1 sub-region; multiple pixels; 
xrt_dem_iterative2.pro, dem_manual.pro; 
with background (subtract has small effect); 
slope increases with AR age 



Table [3] provides a summary of observational measurements in active region cores and 
numerical calculations of the emission measure slope coolward of the peak. The column 'Slopes' 
refers to the range of emission measure slopes coolward of the peak, 'Range of log 10 7" gives 
the temperature range over which the slopes were calculated and 'Reference' the paper in 
which the work is described. The 'Comments' column discusses any particulars that should be 
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considered when making comparisons b etween the differen t studie s. The Markov-Chain Monte 



Carlo (MCMC) method is described by 



Kashyap & Drake 



(I1998|) . xrt_dem _iterative2.pro is a 



SolarSoft (SSW) routine used to generate emission measures constrained at high temperatures 
by Hinode-XRT data, dem_manual.pro is an alternative SSW rout i ne fo r emission measure 



calculations, Pottasch refers to the method described by 



Pottasch 



(11963b and 'Model slopes' 



indicates that the emission measure was calculated from 'perfect' model data. 

Table 4:: Distribution of Observed Emission Measure Slopes 
Coolward of the Peak. 



a < 2.00 2.00 < a < 2.50 2.50 < a < 3.00 3.00 < a < 3.50 a > 3.50 



a 

a - Aa 
a + Aa 



3 
11 



5 
6 



3 
2 
3 



6 
2 
5 



5 
1 

14 



Table |4] shows the distribution of the observed values of a in bins of width 0.5, where we 
have used the average value in cases where the same active region i s observed on the same date. 



Guennou et al. 



(|2012|)have 



The first row assumes that there are no errors in the measurements, 
recently performed a detailed analysis of the effects of realistic measurement errors due photon 
statistics and inaccurate atomic physics. Using a Monte Carlo approach combined with Bayesian 
statistics, they concluded that observationally derived emission measure slopes are uncertain by 
typically +1.0. Atomic physics uncertainties dominate. After consulting several spectroscopists, 
they assumed the following in their synthetic observations: 20% random errors that are different 
for every spectral line (uncertain radiation and excitation rates and atomic structure calculations); 
30% random errors that are the same for every line from a given ion but different for different ions 
(uncertain ionization and recombination rates); 30% random errors that are the same for every ion 
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of a given element but different for different elements (uncertain abundances other than the first 
ionization potential (FIP) bias); and 30% random errors that are the same for all low-FIP elements 
(the uncertain FIP bias). The second and third rows in Table@]indicate the distribution of observed 
a under the assumption that the measurements are in error by Aa = +1 and -1. 



Mulu-Moore et al. I (|201l|) found maximum slopes of 2.3 from their models using coronal 
abundances and suggested that enhanced cooling due to an element population having super- 



coron al abundances in active region cores could give rise to yet steeper slopes. 



Reale & Landi 



(|2012|) demonstrated that new atomic calculations using more sophisticated models and including 
many more emission lines at lower coronal temperatures can yield new estimates of the radiative 
losses that are significantly enhanced, leading to correspondingly steeper emission measure slopes. 
Though the emission measure slopes calculated from forward modeled spectra can be greater than 
the model values in some cases, we note that for low-frequency nano flares it is still difficult to 
obtain slopes that exceed 2.3. There are only three examples in Table[T](Runs [23,40,41]), all for 
relatively low total volumetric heating (< 10 erg cm" 3 ). Nonetheless, forward modeling shows 
that the applicability of low-frequency nanoflares may extend somewhat further into the lower 
end of the range of observational results than may be expected from calculating the slope from 
'perfect' model data alone. 

Our results encapsulate a parameter space within which we investigate variations in the loop 
length, the magnitude and timescale of energy release, and the temporal envelope of the energy 
release mechanism. We now discuss the influence of each parameter in turn, with reference to the 
results listed in Table [T] 

We have chosen lengths such that 20 < 2L < 160 Mm. Assuming approximately 
semi-circular loops then a total length of 160 Mm yields a foot-point separation of about 100 Mm 
or 140", a reasonable upper limit for the diameter of an active region core. The general trend 
that we observe in our numerical results is for the model and forward modeled emission measure 
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slopes to steepen with increasing 2L. This becomes most apparent when comparing particular 
groups of Runs. For example, consider the pairs of Runs [3,10], [4,11], and [5,12] for which 
only the length changes (we discard Run 2 and its corresponding Run 9 because the emission 
measure peaks below 1 MK in Run 2). The slope can be seen to steepen, but of these Runs 
only 11 & 12 may be of interest because the slopes are consistent with the lower end of the 
observed range and ^p ea k ls consistent with the temperature range of peak active region core 
emission (6.45 < log 10 ^p ea k ^ 6.75). Runs 15-26 and 31-39 constitute a more detailed 
study of the relationship between the loop length and the slope. Again, considering equivalent 
triples of Runs where only the length changes (e.g. [15,19,23], [16,20,24], and [31,34,37], etc.) 
we see that the trend is largely maintained, though the slope of the forward modeled emission 
measure for Run 20 is somewhat steep in comparison with Runs 16 & 24. Likewise Run 34. 
These anomalies could likely be corrected by choosing a better estimate for the density in the 
contribution function, bringing the forward modeled slope into closer agreement with the model 
slope. Within these groups of Runs a peak heating rate of 0.5 erg cm" 3 s" 1 yields a T^ Qg ^ that is 
in general too high to be consistent with active region cores for the length range L > 40 Mm. In 
most cases a peak heati ng rate of 0. 1 erg cm " 3 s _1 yields T^ ea ^ in the region of those observed 
(6.40 < log 1() T < 6.65: [Warren et al. 1 1201 21) and so the upper limit to the rate of energy release 
probably lies between 0.1 and 0.5 erg cmr 3 s _1 for timescales of O(100) s (stronger heating may 
be permissible on shorter timescales, we consider this below). Our experiments indicate that, 
regardless of the loop length, low-frequency nanoflares can only account for the lower end of the 
range of observed slopes (< 2.6, assuming no errors) while remaining consistent with the peak 
emission temperature of active region cores. 

We can understand the dependence of the emis sion measure slope on the loop length by 



Bradshaw & Cargill 



(|2010h who showed that the value of the index in 



considering the work of [ 
the scaling law T oc n d depends on the loop length, with smaller values of 6 for longer loops. 
8 itself depends on the ratio of the draining timescale to the radiative cooling timescale in the 
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corona. Longer loops exhibit significant gravitational stratification and lower densities in loop 
apex regions then lead to longer coronal radiative cooling times and consequently smaller 6. The 
peak emission measure arises at the peak density, which occurs approximately at the transition 
from conductive to radiative cooling after the ablative upflows have ceased transporting material 
into the corona. The emission measure slope coolward of the peak is then dominated by strands 
that are cooling by radiation and enthalpy-driven losses. For shorter loops 6 « 2 and the loop 
cools significantly more quickly than it drains. When the strand temperature reaches 1 MK the 
density hasn't fallen too far from its peak. Conversely, for longer loops 6 ~ 1 and the loop cools 
about as quickly as it drains. When the strand temperature reaches 1 MK a greater fraction of the 
material will have been lost through draining than in the short loop case, with a commensurately 
greater change in density. Since EM{T) oc n 2 we can see that a greater density change in long 
loops during the radiative and enthalpy-driven cooling phase will lead to a correspondingly steeper 
emission measure than in short loops in the range 1 MK to Tp &a ^. 

Runs 15-39 were designed to show the effect of varying the peak heating rate on the 
emission measure slope and for this group we have chosen a range of volumetric heating rates 
such that 0.03 < E m < 0.5 erg cm 3 s _1 . The weak trend that emerges is for somewhat steeper 
slopes with increased heating, though it is by no means universal in either the model or forward 
modeled sets of slopes. It is clear from our results that merely changing the heating rate is not 
sufficient to obtain emission measure slopes that are in better agreement with the observed slopes, 
and beyond some upper limit the predicted peak emission measure temperature will exceed the 
temperature of peak emission observed from active region cores. 

The heating timescales that we have chosen for our study lie in the range 10 < t h < 2000 s. 
t h is the key parameter because we are aiming to use the emission measure slope as a diagnostic 
of the heating timescale. Runs 3 - 5, 9 - 12, 23 - 30, and 40 - 45, show the results of varying this 
parameter. Runs 3-5,9-12, and 40 - 45, show a mo ^ e \ generally increasing with t h . There is no 
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clear pattern for ^observed' ^ ut we note mat where particularly steep slopes are found there can 
be considerable scatter in the Pottasch emission measure values and a correspondingly small R 2 
obtained from the linear fit (R 2 is the square of the correlation coefficient obtained from standard 
linear regression). The scatter indicates that the density dependence of the contribution function 
might be an important factor for these emission measures. Runs 23 - 30 (e.g. pairs [23,27], 
[24,28], etc.), for which group the loop length is significantly longer, show an overall decrease in 
^model w * m mcreasm g t h- We see a tendency for observed t0 mcrease with t h (except for the 
weakly heated pair [23,27] indicating a density-related issue) in Runs 23 - 30. 

We can combine the parameters E m and t h into a single parameter: the total volumetric 
energy input / heating rate E H erg cm -3 . There is a clear relationship between 2L, E H and 
log 10 T^eak m Table [T] Longer loops require a weaker total volumetric energy input for 
consistency with t he observed range o f log 10 7"p ea k, as may be inferred from the general coronal 



loop scaling laws (IRosner et al. 



1978|) . There is some evidence in the case of shorter loops to 
suggest that, for Runs with the same (or similar) E H , stronger heating on a shorter timescale yields 
a steeper slope than weaker heating on a longer timescale. For example, consider the pairs [4,6], 
[5,7], and the single Run 8 which has very strong heating on a very short timescale but yields 
a slope at the lower limit of the observed range. The shallow slopes for loops of length 40 Mm 
(Runs 15-18 and 31 - 33) which are heated on the relatively long timescale of 500 s are also 
indicative. There is no evidence of such a relationship for longer loops. Hence, for short loops our 
experiments hint that strong heating may be required for consistency with the observed range of 
a and log 10 7p ea k, which in turn requires a short heating timescale to prevent log 10 Tp ea k from 
overshooting this range. Therefore if one is able to constrain E H and E H0 , from the loop length 
and log 10 7p ea k, then it may also be possible to constrain t h . 

Finally, we can examine the effect of the temporal profile of the energy release on the 
emission measure slope. We choose two profiles: a square / top-hat; and a triangular distribution. 
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To draw out the effects of the shape of the temporal profile we compare pairs of Runs in the 
groups 15 - 26 and 31-39 with the same (or similar) E H . Considering the pairs [16,31], [17,32], 
[20,34], [21,35], [24,37], and [25,38], we find reasonablly good agreement between a mo ^ c \ and 
^observed m eacn comparison and no apparent dependence on the temporal profile. While on 
the one hand it may be a pity that there are no obvious clues to the temporal profile using this 
diagnostic technique, since it is also an unknown in our modeling then it is reassuring to note that 
our predictions of the emission measure slope are not strongly dependent on the form that we 
choose. For a given E H then a square or triangular profile should yield a reasonably consistent 
slope. 

In general, we conclude that any reasonable combination of loop and heating properties can 
explain only the lower end of the range of observed emission measure slope s if th e measurements 



Guennou et al. 



(120121) imply that 



are error free. However, the sizable uncertainties found by 
low-frequency nanoflares might be consistent with a majority of active region cores or with none 
at all. If the loop length and the temperature of peak emission can be measured, and if the slope 
uncertainties can be reduced, then it should be possible to constrain the total volumetric energy 
input E H on each strand. In this case our numerical results show that the emission measure-slope 
diagnostic may be sufficiently sensitive to constrain t# for loops of 2L < 40 Mm. 



4. Summary and conclusions 

We have run a large number of numerical experiments in order to explore an extensive 
parameter space of loop and heating properties so that we may determine what the emission 
measure-slope diagnostic can tell us about the time-dependence of active region core heating. We 
have found that low-frequency nanoflares, where individual magnetic strands cool fully before 
being re-energized, are consistent with only the lower end (1.70 < a < 2.60) of the range of slopes 
reported from observations (1.70 < a < 5.17). These observed slopes have sizable uncertainties, 
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however. a exhibits a fairly strong dependence on the loop length which can be understood in the 
low-frequency nanoflare case because its value is determined by the strands undergoing radiative 
and enth alpy-driven cooling, whe re the relationship between T and n has a strong dependence on 



length (|Bradshaw & Cargill 



|2010|) . In the range consistent with low-frequency nanoflares we have 
found that the emission measure-slope diagnostic may be sensitive enough to yield information 
about the time-dependence of the heating mechanism in the case of short (2L < 40 Mm) loops. 

Considering our results for observed a l° ne > where we have found 0.81 < observed ^ 2.56, 
then Table |4] indicates that approximately 36% of observed active region cores are consistent 
with heating by low-frequenc y nanoflares, assumin g no errors in the slope measurements. The 



uncertainties determined by 



Guennou et al. 



(120 12|) suggest that as many as 77% of cores may be 



consistent with low-frequency nanoflares or as few as none. 

We can place further physical constraints on our results and restrict the parameter space by 
considering the nature of the heating mechanism. If we assume that the heating arises from the 
impulsive release of energy from the magnetic field then E H must be limited to the amount of free 
energy available in the field. The free magnetic energy density is given by: 



(i) 



where B p is the potential component of the magnetic field and e is a parameter that indicates the 
level of stress suc h that B x = eB p is the stress component. A typical value is e = 0.3, with an 



upper limit of 0.5 (|Dahlburg et al. 



2005). By setting E B = E H we can determine the magnetic 



field strength B needed to provide the total volumetric heating, where B 2 = B 2 + B 2 . We must bear 
in mind that not all of the free energy may be released from the field at one time, which means 
that the magnetic field strength estimated from Equation Q] may be considered a lower limit. 



Mandrini et al. 



(|2000l) studied how the average values of B and B 2 along a field line depend 
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on the field line length, 2L. Table [5] gives the results for several observed active regions based 
on equation 9 of their paper. The active regions (numbers 1-6 in the paper) cover a range of age 
and complexity. Only field lines with photospheric field strengths between 100 and 500 G were 
used to derive the results. Table[5]lists < B 2 > l/2 for the four loop lengths used in our nanoflare 
simulations. Values are given for the weakest and strongest active region and the average of the 6 

Table 5:: Average magnetic field strength (<B 2 > 112 ^ versus 
loop length. 



2L 


^min 


#max 


#avg 


[Mm] 


[G] 


[G] 


[G] 


20 


127 


211 


167 


40 


83 


189 


136 


80 


42 


150 


94 


160 


18 


89 


51 



Table [6] shows a subset of the experiments from our study that most closely match 
observations and also satisfy the magnetic energy constraint. The criteria for selection are that 
6.45 < log 10 7"p ea k ^ 6.75, ^observed - ^ an ^ B < Bavg, where 5 a vg is the length-dependent 
value in the last column of Table \5\ Regular font indicates that the most restrictive energy 
constraint e = 0.3 satisfies these criteria, and italics indicate that only the less restrictive constraint 
e = 0.5 can satisfy these c riteria. The Runs are listed in order of increasing Q'observed- Of 



the original 14 Runs from 



Bradshaw & Klimchuk 



(1201 lh . only [8,11,12,13,14] are eligible 



for inclusion in Tabled We can see that the maximum value of E H must be somewhere in 
the region of 50 erg cm" 3 to avoid unphysical values of B, though exceptions are possible 
(e.g., Run 18), especially in magnetically strong active regions (third column of Table [5]). In 
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the case of short loops (< 20 Mm) this limit indicates strong volumetric heating (to reach 
6.45 < log 10 7"p ea k < 6.75) and short timescales (to satisfy E H < 50 erg cm 3 ). Run 8 provides 
a good example. For longer loops weaker volumetric heating over extended timescales, in 
comparison with shorter loops, yields emission measure peaks, observed slopes and magnetic 
field strengths in the right range. 

Table 6:: The Subset of Numerical Experiments that Most 
Closely Match Observations. 



Run # 2L E H a mo del ^observed B ( 6 = °- 3 ) B ( e = °- 5 ) 





[Mm] 


[erg cm 3 ] 






[G] 


[G] 


8 


20 


50 


1.79 


1.73 


123 


79 


37 


160 


15 


1.98 


1.76 


68 


43 


24 


160 


15 


2.01 


1.77 


68 


43 


14 


80 


30 


2.01 


1.80 


96 


61 


36 


80 


50 


1.81 


1.82 


123 


79 


12 


80 


30 


1.90 


1.85 


96 


61 


19 


80 


7.5 


1.40 


1.85 


48 


31 


44 


80 


30 


2.05 


1.90 


96 


61 


18 


40 


125 


1.95 


1.93 


195 


125 


42 


80 


6.0 


1.91 


1.93 


43 


27 


13 


80 


10 


1.90 


2.12 


55 


35 


20 


80 


15 


1.53 


2.16 


68 


43 


34 


80 


15 


1.81 


2.16 


68 


43 


45 


80 


60 


2.04 


2.16 


135 


87 


11 


80 


10 


2.25 


2.22 


55 


35 


43 


80 


12 


1.95 


2.28 


60 


39 
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Table 6:: The Subset of Numerical Experiments that Most 
Closely Match Observations. 



Run# 


2L 


E H 


a model 


•^observed 


B(e = 0.3) 


B(e = 0.5) 




[Mm] 


[erg cm 3 ] 






[G] 


[G] 


23 


160 


7.5 


1.68 


2.47 


48 


31 


41 


80 


3.0 


1.46 


2.56 


30 


19 



In summary; we can use the loop length and the magnetic field strength to constrain E H and, 
of the components of E H , we can use 7"p ea k to constrain E H0 and then E H /E H o yields t H - Our 
studies show that the emission measure slope is essentially independent of the temporal envelope 
chosen for the heating and so t h is the key parameter. For low-frequency nanoflares we find 
predicted slopes in the range a < 2.6. 

Finally, we can place limits on the slope that can be obtained in the case of heating by 
low-frequency nanoflares. The emission measure slope coolward of the peak is determined by the 
strands undergoing radiative cooling and enthalpy-driven draining. The emission measure in a 
particular temperature range depends upon two factors: (a) the density of the strands as they cool 
through that range; and (b) the number of strands in that range. A characteristic emission measure 
can be calculated by weighting the density-dependent emission measure by the number of strands. 
Since the number of strands in a particular temperature range depends upon how quickly the 
plasma cools through that range, At(T), we can write 



(EM(T)) = — !— f EM(T)At(T), (2) 
T cool J 



where t coo j is the total cooling time. We know that 
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EM(T) = n 2 D, 



(3) 



where D is the line-of-sight depth of the emitting plasma, and 



Ar(D 



n 2 A(T) 



where A(T) = xT 



(4) 



To address the first of the above factors, (a), we r equire an expression for the density in terms of 



the temperature of the cooling strands. Following 



Bradshaw & Cargill 



(2010) we can write 



T oc n 5 where 5 = (y - 1) + 

r CR 



(5) 



where y = 5/3, Ty is the coronal draining timescale and t cr is the coronal radiative loss timescale. 
We then have 



n oc T 6 , (6) 

and since EM(T) oc n 2 we can write 

EM(T)ocT Z *. (7) 

To address the second factor, (b), we require an expression for the cooling timescale in terms of 
the temperature. Using Equation @] we can see that 

j\-b 

At(T) oc (8) 

n 



and substituting Equation [6] gives 
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At(T) oc T l - b ~. (9) 



We can see from Equations [2l|7] and |9] that 



(EM(T)) ocTh l - h - L t 

and so 



(EM{T)) oc T a where a = - + 1 - b. (10) 

6 

The slope of the emission measure, a, in any temperature range then depends upon the parameter 
8 and the slope of the radiative loss function b. 

We can place an absolute upper limit on a by recalling that in the limit of pure enthalpy-driven 
cooling (jv — > and negligible corona radiative losses) 6 = y — 1 = 2/3 and 



amax = ^ ~ b - 



(ID 



Typical values of b in t he temperature range over w hich emissio n measure slopes are ca lculated 



are given by b = -1/2 (|Bradshaw & Cargill 



2010) or b = -3/2 (Klimchuk et al 



Bradshaw & Cargill 



2008). These 



d20JD) showed 



yield a(b = -l/2) max = 3.0 and a(b = -3/2) ma x = 4.0. 
that 6 decreases with increasing loop length and to approach 6 = 2/3 would require an unfeasibly 
long loop for our study. We are limited to loops of app roximately 160 Mm i n tota l length to 



Bradshaw & Cargill 



(2010) shows that 



be consistent with active region core scales. Table 2 of [ 
6(L = 160 Mm) » 1.50, which yields a(b = -1/2) = 2.17 and a(b = -3/2) = 3.12. In the 



case of shorter loops Table 2 of 



Bradshaw & Cargill 



(12010 ) shows that 6 ~ 2.00, which yields 
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a(b = -1/2) = 2.00 and a{b = -3/2) = 3.00. We note that lSturrock et al. I (119901) performed a 
related analysis and found that EM(T) oc T l ~ b . 

In consequence, we conclude that it is not possible for the low-frequency nanoflare scenario 
to yield emission measure slopes coolward of the peak that exceed these values of a, depending 
upon the radiative loss function. In the case of b = -l/2(-3/2), approaching a = 3(4) would 
require loops of unphysical length for active region cores. This explains why our experiments 
are limited to slopes of about 2.6, because the maximum length in our study is 160 Mm 
(2.00 < a < 3.00 for the values of b considered). Therefore, we determine that low-frequency 
nano flares cannot explain emission measure slopes in active region cores where a > 3.00. 

A partial solution to the difficulty of obtaining sufficiently steep emission measure slopes 
with low-frequency nanoflares might be found from the super-position of coronal loops. Consider 
a scenario in which a number of particularly long loops (2L > 160 Mm, for example) cross 
the active region core along the observer's line-of-sight, but have foot-points located outside 
the core. The contribution to the emission measure of such loops might serve to raise the limit 



on a due to their smaller values of 6, though they are stil 
which yields atnax less t han the steepest slopes found by 



Schmelz & Pathak 



bound by the lower limit 6 = 2/3 



Warren et al. 



(120121) (a = 4.50) and 



(|2012|) (a = 5.17). One might argue that the higher densities of short loops 
(shorter than a gravitational scale height) means they would make the dominant contribution to 
the line-of-sight emission and hence to a, but their radiative lifetimes are also short meaning that 
fewer strands contribute to the emission in a given temperature range (equation |2]). Consequently, 
we might expect long loops and short loops to have similar emission measures. We can resolve 
the question of which loops make the dominant contribution by noting that density is strongly 
dependent on length, and length only enters the weighted emission measure via the parameter 6. 
We have already demonstrated that ornax is smaller for shorter loops and restricted to an upper 
limit, that lies below many observationally measured values, for longer loops. 



-32- 



In practical terms, then, it seems that to find slopes greater than 2.6 for current radiative loss 
calculations we must appeal to some other heating scenario. Low-frequency nanoflares are one 
limit, where the time between heating events on a single strand (the re-energization timescale) is 
longer than the cooling / draining timescale. The opposite limit is steady heating, where the strand 
has no time to cool before the onset of the next heating event. In order to steepen the slope it is 
necessary to increase the amount of hot material in the region of ^p ea k relative to the amount of 
material near 1 MK. In the next paper of this series we will investigate repeating nanoflares on a 
single strand, where the re-energization timescale is less than the cooling / draining timescale, in 
order to determine whether this mechanism can maintain enough hot material, relative to cooler 
material, to explain the steepest observed slopes (e.g. a > 3). 
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